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With A Multi-Phase Transport (AMPT) model we investigate the relation between the magnitude, 
fluctuations and correlations of the initial state spatial anisotropy e n and the final state anisotropic 
flow coefficients v n in Au+Au collisions at ^s NN = 200 GeV. It is found that the relative eccentricity 
fluctuations in AMPT account for the observed elliptic flow fluctuations, both are in agreement with 
the elliptic flow fluctuation measurements from the STAR collaboration. In addition, the studies 
based on 2- and multi-particle correlations and event-by-event distributions of the anisotropies 
suggest that the Elliptic-Power function is a promising candidate of the underlying probability 
density function of the event-by-event distributions of e n as well as v n . Furthermore, the correlations 
between different order symmetry planes and harmonics in the initial coordinate space and final 
state momentum space are presented. Non-zero values of these correlations have been observed. 
The comparison between our calculations and data will, in the future, shed new insight into the 
nature of the fluctuations of the Quark-Gluon Plasma produced in heavy ion collisions. 

PACS numbers: 25.75.Ld, 25.75.Gz 


I. INTRODUCTION 

One of the fundamental questions in the phenomenol¬ 
ogy of Quantum Chromo Dynamics (QCD) is what the 
properties of matter are at the extreme densities and tem¬ 
peratures where quarks and gluons are in a new state of 
matter, the so-called Quark Gluon Plasma (QGP). Colli¬ 
sions of high-energy heavy-ions, at the Brookhaven Rela¬ 
tivistic Heavy Ion Collider (RHIC) and the CERN Large 
Hadron Collider (LHC), allow us to create and study the 
properties of such a system in the laboratory. 

The azimuthal anisotropy in particle production is, at 
these energies, an observable which provides experimen¬ 
tal information on the Equation of State and the trans¬ 
port properties of the QGP. This anisotropy is usually 
characterized by the Fourier flow-coefficients [1], 

v n = (cos[n((£> — x L ra )), or equivalently 

v n = (e inv e~ in ^ n ), ( 1 ) 

where is the azimuthal angle of the particles, ’P Tl is 
the n th -order flow plane (or named final state symmetry 
plane) angle and ( ) denotes an average over the selected 
particles and events. 

In the last decade, the elliptic flow i >2 [2-8], which is 
considered to correspond to the elliptical shape of the 
spatial overlap region in the system created in the colli¬ 
sions [9], has received a lot of experimental and theoret¬ 
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ical attention. For a recent summary see [10-15]. More 
recently, higher odd and even anisotropic flow coefficients 
are found to be also very important [16]. Hydrodynamic 
calculations predict that these higher harmonics, such as 
the triangular flow U 3 , are more sensitive to the shear 
viscosity to entropy density ratio rj/s of the QGP than 
V 2 [17]. Furthermore, it is realized that the correlations 
between the symmetry planes and flow harmonics are 
sensitive to both the initial state and r\/s [18-20]. The 
combined analysis of both e n and v n distributions for a 
single harmonic, and the correlations between different 
orders of symmetry planes (direction of the flow vector) 
and flow harmonics (magnitude of the flow vector) could 
yield powerful constraints on both initial conditions and 
properties of the QGP. However, so far the investigation 
of the above mentioned combined analysis using a trans¬ 
port model was still lacking. 


In this paper, we present the calculations of initial 
state anisotropies and final state anisotropic flow in 
Au+Au collisions at ^/s NN = 200 GeV using the AMPT 
model [2 ]. We also investigate the relation between the 
magnitude, fluctuations and correlations of the initial 
state spatial eccentricity and final state anisotropic flow 
coefficients. In addition, the correlations between differ¬ 
ent order symmetry planes and harmonics will be inves¬ 
tigated in both initial and final state, which will help us 
understand how they are modified during the expansion 
of the system. 
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II. A MULTI-PHASE TRANSPORT MODEL event weight proposed in [32]: 


A Multi-Phase Transport (AMPT) model [2 ] with a 
so-called string melting scenario has been used for these 
studies. The model consists of four main stages: ini¬ 
tial conditions, partonic interactions, hadronization and 
finally hadronic rescattering. 

The initial conditions, which include the spatial 
and momentum distributions of minijet partons and 
soft string excitations, are obtained from the HIJING 
model [22]. The strings are converted into partons and 
the next stage, which models the interactions between all 
the partons, is based on ZPC [23]. The partonic cascade 
model ZPC presently includes only 2-body processes 
with cross sections obtained from pQCD with screen¬ 
ing masses. In ZPC, the default value of the cross sec¬ 
tion is 3 mb. The transition from partonic- to hadronic- 
matter is modeled by a simple coalescence model, which 
combines two quarks into mesons and three quarks into 
baryons [24]. Finally, to describe the dynamics of the 
subsequent hadronic stage, a hadronic cascade based on 
the ART model [25] is used. In this analysis, we used the 
default input parameters of AMPT with string melting 
suggested in [21]. The results are presented as a function 
of centrality, which determined by impact parameter b, 
as used in [26, 27]. The possible effects of the fluctuations 
of impact parameters will be discussed in next section. 


III. ANALYSIS METHOD AND DEFINITIONS 

In this paper the anisotropic flow is calculated us¬ 
ing the 2- and multi-particle cumulants method [28, 29], 
which was widely used at RHIC [30] and at the LHC [ 6 , 
31]. In this method, both 2- and multi-particle azimuthal 
correlations are analytically expressed in terms of a Q- 
vector, which is defined as: 


M 

Q n = 

i= 1 


( 2 ) 


where M is the multiplicity of the selected particles and 
ip is their azimuthal angle. 

The single-event average 2-particle azimuthal correla¬ 
tions can be calculated via: 


_ \Q n \ 2 -M 
' ; M(M - 1)' 


(3) 


From this the event averaged 2-particle correlations, and 
the 2 -particle cumulants can be obtained using: 


Cn{2} = ((2)) 


X) events (^( 2))i(2)i 

X ]events (^^{2) )i 


(4) 


where W^ 2 ) is the event weight. To minimize the effect of 
the varying multiplicity in certain centrality class deter¬ 
mined by b, we use the number of particle pairs as the 


W {2 ) = M{M - 1 ) (5) 

The anisotropic flow from 2-particle cumulants, denoted 
as v n {2} 1 is finally obtained from: 

v n {2} = \J c n {2}. (6) 


Unfortunately, the v n {2} contains contributions from 
so-called non-flow effects, which are additional azimuthal 
correlations not associated with the common symmetry 
planes, e.g. resonance decays, jet fragmentation, and 
Bose-Einstein correlations. They can be suppressed by 
appropriate kinematic cuts and therefore one can intro¬ 
duce a gap in pseudorapidity between the particles used 
in the 2-particle Q-cumulant method [33]. For this we 
divide the whole event into two sub-events, A and B , 
which are separated by a pseudorapidity gap \Ar]\. This 
modifies Eq. (3) to: 


(2) A r] 


Qn ■ QZ* 

M a ■ M b ’ 


(7) 


where Q£ and Q ® are the flow vectors from sub-event A 
and B 1 with M A and Mb the corresponding multiplici¬ 
ties. The event weight from Eq. (5) in this case becomes: 


W {2)av = M a -M b . (8) 

Finally, inserting Eqs. (7) and ( 8 ) into Eq. (4), the v„ 
from a 2 -particle cumulant with a A 77 gap is given by: 

Vn{2,\A v \} = ^/((2)) Ari . (9) 

Instead of using kinematic cuts, the collective nature 
of anisotropic flow itself can be exploited to suppress 
non-flow contributions. Using multi-particle instead of 
2 -particle cumulants the aforementioned non-flow effects 
are strongly suppressed and no additional kinematic cuts 
are required. The v n calculated using cumulants are de¬ 
noted as v n {k}, were k is 2 ,4, 6 ..., m for the m-particle 
cumulant. Following the Q-cumulant method [29], the 
single event average 4- and 6 -particle correlations can be 
calculated as: 

(4) = [ IQnl 4 + \Q 2 n\ 2 - 2 • Re (i Q2nQnQ*n ) 

— 2[2(M — 2) • \Q n \ 2 — M(M — 3)] ] 

/[M(M - 1 ){M - 2){M - 3)], 

(6) = [ IQnl 6 + 9\Q2n\ 2 \Qn\ 2 - 6 ■ Re ( Q2nQnQ* n QnQn) 

+ 4 • Re (Q 3n Q* n Q* n Q* n ) - 12 • Re (Q 3n Q* 2n Q* n ) 

+ 18(M - 4) • Re ( Q2nQ* n Qn ) + 4|Q 3 n| 2 

— 9(M - 4)(\Qi\ + \Q'2n | 2 ) + 18(M - 2 ){M - 5)|Q„| 2 

— 6M (M — 4) (M — 5) ] 

/ [M(M - 1 )(M - 2)(M -3)(M - 4 ){M - 5)] 

( 10 ) 
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Then the multi-particle cumulants are obtained from: 

Cn{4} = ((4»-2((2)) 2 , 

Cn{6 } = ((6))-9((2)).((4)) + 12((2)) 3 , 


here (( )) denotes the average over all particles over all 
events. In the end the v n from 4- and 6 -particle cumulant, 
denoted as rr n {4} and u n { 6 }, are obtained using: 

«n{4} = \/~c n { 4}, 
v n {6} = \/ c ra { 6 }. 

The 2- and multi-particle cumulants have different con¬ 
tributions from flow fluctuations. The contribution from 
flow fluctuations is positive for the 2 -particle cumulant 
and negative for the multi-particle cumulant [34, 35]. 
When non-flow effects are negligible for the 2-particle 
cumulant, and if a Vn <C v n , the v n from the cumulants 
are up to order cr 2 ^ given by: 

v n {2} 2 « v 2 n +al n , 

v n {4 } 2 «<-<, (13) 

m 2 —2 2 

~ V n - (T Vn , 

where v n and a Vn are the mean and standard deviation 
of the v n distributions, respectively. In the special case 
that the underlying probability density function ( p.d.f .) 
of v n is described by a Bessel-Gaussian function, Eq. (13) 
is an exact solution, independent of the magnitude of the 
flow fluctuations. This shows that for a Bessel-Gaussian 
p.d.f. all the multi-particle cumulants will be identical. 

It is thought that the development of anisotropic flow 
is controlled by the anisotropies in the pressure gradients 
which in turn depend on the shape and structure of the 
initial density profile. The latter can be characterized, in 
analogy with the flow Fourier coefficients and flow angles 
of Eq. (1), by a set of harmonic anisotropy coefficients e n 
and associated symmetry angles <f>„ of the initial spatial 
distribution: 


AMPT. Figure la also shows the calculated r; 2 {4} and 
•u 2 {6}. Because the non-flow contributions to the inte¬ 
grated flow are very small, the observed significant dif¬ 
ference between u 2 {2} and u 2 {4} is due to elliptic flow 
fluctuations. In addition we observe that there is very 
good agreement between t> 2 {4} and rr 2 { 6 }. This seems to 
agree with the expectation that the underlying p.d.f. is a 
Bessel-Gaussian function, which predicts u 2 {4} = t; 2 { 6 }. 
However, it is still important to point it out that in fact 
i> 2 {4} ~ u 2 { 6 } is valid irrespectively of the details of the 
underlying model of flow fluctuations, under the assump¬ 
tion <x„ n -C uf [ ]. 


>C " 0 . 08 b 


☆ v 2 {2} 

A v 2 {2,|Ar||>1} 
f> v 2 {est} 

■ v 2 {4} 

♦ v 2 {6} ^ 



30 40 50 60 

centrality percentile 




inG> n 


f r dr df>r n e m ^ e(r,<f>) 

J r dr d4>r n e(r,(j)) ™ 


(14) 


where e(r, f) is the initial energy density distribution in 
the plane transverse to the beam direction. Also these 
anisotropy coefficients can be calculated using cumulants 
and the same relation between e n and a Eri follows as in 
Eq. 13, replacing v n with e„. 


IV. RESULTS AND DISCUSSION 


FIG. 1: (Color online) Centrality dependence of (a) V 2 and 
(b) £2 in Au+Au collisions at ^s NN = 200 GeV in AMPT. 
For the definition of the symbols see text. 


From u 2 {2, \Ar/\ > 1} and v 2 {4}, we can estimate u 2 
in Eq. (13) using: 


w 2 {est} = 


u 2 {2, |Ary| > l } 2 +u 2 {4 } 2 


(15) 


The centrality dependence of u 2 from the AMPT model 
calculations is shown in Fig. la. For the 2-particle cu¬ 
mulant calculations, we plot both u 2 {2} from Eq. ( 6 ) 
and r; 2 {2, |Ary| > 1} from Eq. (9). It is seen that the 
u 2 {2, | A?y| > 1} is compatible with u 2 {2}, which indicates 
that short range non-flow contributions are very small in 


Figure lb shows e 2 as function of centrality. Here e 2 
can be calculated from the initial spatial parton distri¬ 
butions in the AMPT model. The relative fluctuation of 
£ 2 , named cr e2 , is shown by the green band as the un¬ 
certainty of e 2 . The £ 2 ( 2 } and £ 2 {4} correspond to the 
2- and 4-particle cumulant definition but are evaluated 
using the initial spatial coordinates. We see that £ 2 {2}, 
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e 2 and £2 {4} increase monotonically up to 60% centrality 
percentile. In contrast, the n 2 , plotted in Fig. la, starts 
to saturate or decrease from 40% centrality percentile. 
This difference between the centrality dependence of v 2 
and £ 2 is generally understood to be due to the fact that 
the efficiency of converting the initial eccentricity into 
final elliptic flow decreases towards peripheral collisions 
because of the decreasing number of reinteractions in a 
smaller system. 



FIG. 2: (Color online) Centrality dependence of (a) ellip¬ 
tic flow fluctuation and (b) relative elliptic flow (eccentricity) 
fluctuations in Au+Au collisions at 200 GeV in AMPT String- 
Melting. (a) The full stars are the estimated elliptic flow fluc¬ 
tuation <t„ 2 and the dashed line is the relative eccentricity 
fluctuations scaled by a factor of 0.15; (b) the estimated rela¬ 
tive elliptic flow fluctuation R V2 4 , the relative eccentricity and 
true eccentricity fluctuation for the AMPT model are shown 
using full diamonds, a dashed line and solid line. In addition, 
the measurement from the STAR Collaboration are plotted 
using full triangles and the shadow shows its uncertainty [30]. 


From the calculated anisotropic flow, we can investi¬ 
gate the relation between the fluctuations in e n and v n . 
In Fig. 2a we plot the centrality dependence of the elliptic 
flow fluctuations a V2 estimated via 


v 2 { 2 , \At)\ > l} 2 - u 2 {4} 2 


(16) 


We find that a V2 increases toward peripheral collisions 


and its magnitude is significant compared to the magni¬ 
tude of v 2 . Assuming v 2 oc £2 and using Eq. (13), we 
also have v 2 {2} oc £ 2 {2} and a V2 oc cr £2 . Hence, fig¬ 
ure 2a also shows the scaled cr E2 calculated as a e2 = 
\J (e 2 {2} 2 — £ 2 {4} 2 )/2. It shows that the eccentricity 
fluctuations describe the centrality dependence of the el¬ 
liptic flow fluctuations quite well, which might indicates 
that they are the dominant contribution for the observed 
elliptic flow fluctuations. 

The estimated relative flow fluctuations, can be calcu¬ 
lated using: 


Ry = 

2,4 


/n 2 { 2 , | A 77 I > l } 2 - v 2 { 4} 2 
v 2 { 2 , \Ar]\ > l} 2 + v 2 { 4} 2 


— , (17) 

V2 


and its centrality dependence is shown in Fig. 2b. The 
magnitude of the relative fluctuations ranges from 0.4 to 
0.6 in central to mid-peripheral collisions. This shows 
that the assumption a V2 <C v 2 does not hold for these 
centralities. 

The estimated relative eccentricity fluctuations R e 
can be calculated analogously to the R v : 


Re 


2,4 


£ 2 {2} 2 - £ 2 {4} 2 _ Oe2 
£ 2 {2} 2 +£ 2 {4} 2 £ 2 ’ 


(18) 


and compared to R Vo . This comparison does not depend 
on the assumption that the relative fluctuations are small 
or the underlying p.d.f. is a Bessel-Gaussian. We notice 
that the i? £a (also plotted in Fig. 2b) is compatible with 
R v for central and mid-central collisions. To test if this 

^ 2,4 

is only accidental we also study the v 2 in AMPT with 
an approximately three times larger partonic cross sec¬ 
tion (10 mb). The magnitudes of both v 2 {2} and v 2 {4} 
increase significantly, however, the consistency between 
R e2 4 and R V2 4 continues to hold. This is expected if 
the relative elliptic flow fluctuations depend only on the 
eccentricity fluctuations, which again shows that eccen¬ 
tricity fluctuations play an important role in the devel¬ 
opment of elliptic flow fluctuations in the AMPT model. 

In addition, we calculate directly from the initial state 
of the AMPT model the true relative eccentricity fluc¬ 
tuations <j e2 /e 2 (plotted as the the green solid line in 
Fig. 2b). The results are consistent with the estimated 
relative eccentricity fluctuations in the 20-50% centrality 
percentile, while they deviate for central and peripheral 
collisions. Because the assumption a e2 <C £2 is not sat¬ 
isfied over the whole centrality range this indicates that 
for 20-50% the fluctuations might be approximately de¬ 
scribed by a Bessel-Gaussian but to describe the overall 
event-by-event v 2 distributions, we need to search for a 
better candidate of the underlying p.d.f. To see how these 
fluctuations compare to experimental data we compare to 
STAR measurements [30]. The relative elliptic flow fluc¬ 
tuations measured in Au+Au collisions at 200 GeV are 
plotted in Fig. 2b and are in very good agreement with 
the AMPT model calculations. 
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Hydrodynamic calculations have shown that in a given 
event, the v n is a linear response to the initial anisotropy 
( v n = k n e n ), for n = 2, 3 [36]. In the discussion above, 
we found that the elliptic flow and its fluctuations can be 
nicely described by the initial eccentricity together with 
its fluctuations. In order to further constrain the under¬ 
lying p.d.f. of the v n and e n distributions, we study the 
event-by-event distributions of £„, simulated by AMPT. 

In this study, the event-by-event e n distributions are 
investigated in the selected centrality classes. To bet¬ 
ter extract the information on the underlying p.d.f. of 
the e n distributions, several candidate functions are used 
in this paper. One popular parameterization of the e n 
distribution is the Bessel-Gaussian distributions [37]: 

p< £ ») = <“) 

where £o is the anisotropy w.r.t. the reaction plane and 
a is the fluctuation in the spatial anisotropy. It was al¬ 
ready shown in this paper as well as in previous stud¬ 
ies [37, 38] that a Bessel-Gaussian distribution nicely 
describes the £2 distributions for mid-central collisions. 
However, it is not expected to work perfectly in periph¬ 
eral collisions, due to the lack of a constraint that £2 < 
1 in each event [39]. To fix this problem, a simple one- 
parameter power-law distribution [40] 

p(£„) = 2 ae n (1 - 4 ) Q_1 , ( 20 ) 

was proposed to parameterized the fluctuations driven 
anisotropies. Here a quantifies the fluctuations. 

Recently, a new function, named “Elliptic Power” dis¬ 
tribution was proposed in [39] as: 


P(£n) 


OL Sfi 
7T 


(1 


- £n 




f 2 * (1-Q al d0 
Jo (1 - £0 Sn COS 4>) 2a+1 ’ 


( 21 ) 

where a and £q have the same meaning as above. In our 
paper, these three candidates of the underlying p.d.f .s 
are fitted to the £ n distributions in AMPT model. 

Figure 3 shows the £2 distributions in the centrality 
range 0-5%, 30-40% and 60-70% from the AMPT initial 
state. We fit the three distributions with the Elliptic- 
Power (solid lines), Power-Law (dash line), and Bessel- 
Gaussian functions (dot-dash line). It is shown that in 0- 
5% these three functions give consistent results and they 
all fit the £2 distributions quite well. It is understood that 
in the case where £0 <C 1 (small anisotropy from the re¬ 
action plane) and where a > 1 (strong fluctuations), the 
Elliptic Power distribution turns into a Bessel-Gaussian 
distribution; while with £0 = 0 (anisotropy is solely due 
to fluctuations), the Elliptic Power distribution reduces 
to a Power-Law distribution. The nice agreement be¬ 
tween the three functions shows that the £ 0 must be 
very small, which means that the eccentricity is gener¬ 
ated mainly by fluctuations. In 30-40%, it is clear that 
only the Elliptic-Power function agrees with the £2 dis¬ 
tributions. For more peripheral collisions, the Elliptic- 
Power function describes the £2 distributions still quite 



FIG. 3: (Color online) £2 distributions in the AMPT initial 
state 


well, while the Bessel-Gaussian can’t quantitatively re¬ 
produce the £2 distributions. The Power-Law function, 
on the other hand, is not expected to describe the £2 
distributions in non-central collisions, since the eccen¬ 
tricity from the reaction plane (e 0 ) is missing and for 
non-central collisions £q is non-zero and important to de¬ 
scribe the £2 distributions. Furthermore, we calculate 
a E2 /e 2 from the Elliptic-Power function with parameters 
extracted from the fits of the £2 distributions. The result 
shown in Fig. 2b is in nice agreement with the er E2 /£2 
directly from the AMPT initial state. For describing all 
the centralities the Elliptic-Power function gives the best 
description of the £2 distributions. Unfortunately, the ex¬ 
pectation of a 1% difference between the 4- and 6 -particle 
cumulants [39] from the Elliptic-Power function can not 
be tested with the current statistics. 

It is shown in hydrodynamic calculations that the 
higher harmonic flow coefficients are more sensitive to 
both the kinematic viscosity and to the initial geometry 
and its fluctuations [17]. In absence of fluctuations u 3 is 
zero due to symmetry constraints. Figure 4a shows the 
centrality dependence of u 3 obtained for AMPT events 
from different analysis methods. If V 3 is due to event- 
by-event fluctuations in the initial spatial density distri¬ 
butions we would expect that there is no, or only a very 
small, correlation between both the reaction plane angle 
Wrp (spanned by the impact parameter vector and the 
beam direction) and the angle of the 2 rad -order flow plane 
'Fo with respect to the 3 rd -order flow plane H/ 3 . The cor¬ 
relations between 'll 3 and Trp can be studied via: 

«3{^RP} = ((COS 3(<p - Tr P ))) 

= ((cos3(^ - \F 3 ) cos 3(^ 3 - T rp ))) (22) 
= (v 3 (cos3(4/ 3 - ’Frp))). 

In Fig. 4a we observe that u 3 {\E , rp} is consistent with 
zero. It indeed shows that there is no correlation (or 
an extremely weak correlation) between \l / 3 and H/rp in 
the presented centrality range. We also see that both 
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FIG. 4: (Color online) (a) Centrality dependence of u 3 in 
Au+Au collisions at 200 GeV from AMPT model calculations. 
The open stars (red), triangles (black), circles (gray), dia¬ 
monds (blue) and crosses (black) for v 3 {2}, r> 3 {2, |A^| > 1}, 
vfl'IC} and u 3 {'Frp} respectively, (b) The corresponding £3 
as a function of centrality. 


the v 3 {2} and fs{4} are nonzero and show only a weak 
centrality dependence, which is in qualitative agreement 
with earlier observations based on hydrodynamic model 
calculations [17]. 

In Fig. 4b the corresponding £r 3 {2}, £ 3 {4} and £3 are 
plotted which, in contrast with ^ 3 {2} and u 3 {4}, in¬ 
crease with increasing centrality by a factor two, this 
difference might be due to strong viscous damping ef¬ 
fects on v 3 compared to v-i [17]. The expectation that 
£ 3 { 4 } = £ 3 { 4 'rp} = 0 if the p.d.f. would be a Bessel- 
Gaussian function is not observed in the 5-40 % central¬ 
ity percentile range, as is shown in Fig. 4 [5 ]. Similar 
non-zero £ 3 { 4 } were found for wider centrality ranges 
from both MC-Glauber and MC-KLN calculations [41]. 
In fact, non-zero values of not only r; 3 {4} but also U 3 { 6 } 
have been measured in experiments [42]. Compared to 
the Bessel-Gaussian distribution, which was the widely 
used description of the underlying p.d.f, the Elliptic- 
Power function seems to gives an improved description of 
tr E3 /£3 and agrees with a non-zero value of multi-particle 
cumulants of £3 and V3. 

For completeness, Fig. 6 shows the £3 distributions in 



FIG. 5: (Color online) a e3 /£3 distributions in AMPT. 



FIG. 6: (Color online) £3 distributions in the AMPT initial 
state 


the centrality range 0-5%, 30-40% and 60-70% from the 
AMPT initial state. Similar to the study of the £2 distri¬ 
butions, we fit all the £3 distributions with the Elliptic- 
Power (solid lines), Power-Law (dash line), and Bessel- 
Gaussian functions (dot-dash line). It is seen that the 
Bessel-Gaussian function reproduces the £3 fairly well 
except for peripheral collisions. Nevertheless, the expec¬ 
tation of a Bessel-Gaussian, which is £ 3 ( 4 } = £ 3 ( 6 } = 
£ 3 {BP} = 0 (as well as r> 3 { 4 } = r; 3 { 6 } = ^{RP} = 0 ), 
does not agree with the non-zero £ 3 ( 4 } and ^{4} pre¬ 
sented above. The Bessel-Gaussian function is therefore 
not a candidate of the underlying p.d.f. of v n and e n . 
On the other hand, since triangularity is expected to be 
solely created by initial geometry fluctuations, its distri¬ 
butions should be well reproduced by a single-parameter 
Power-Law function [43]. Indeed nice agreements be¬ 
tween the fits of the Power-Law and the £3 distributions 
are observed for the presented centralities. When test¬ 
ing the two parameter Elliptic-Power function, it turns 
out, as expected from the nice fit of the Power-Law func¬ 
tion, that the parameter £q (triangularity w.r.t. reaction 
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plane) is very close to 0. The nice descriptions of the £3 
distributions by the Power-Law function and the Elliptic- 
Power function with parameter e 0 ~ 0 confirms that the 
contributions of the reaction plane to £3 is very weak (or 
zero). This also agrees with our results of u 3 {RP} = 0, 
displayed in Fig. 4, and with the experimental measure¬ 
ments of ALICE [31]. Furthermore, we notice that the 
true a e3 /£3 from the AMPT initial state shows a decreas¬ 
ing trend from central to peripheral collisions. This de¬ 
creasing trend is captured by the Elliptic-Power function 
quite well, and clearly disagrees with the expectation of 
a Bessel- Gaussian type p.d.f., which predicts a constant 
value of a/ 4 / 7 t — 1 for the entire centrality range [10, 38]. 
It is shown that the Elliptic-Power function gives a bet¬ 
ter description of the £2 and £3 distributions simultane¬ 
ously, matches the multi-particle cumulants observables 
and the correlations w.r.t. the reaction plane, and also 
matches the experimental measurements [31]. This shows 
that the Elliptic-Power distributions could be a promis¬ 
ing candidate of the underlying p.d.f. of v n (and e n ), 
which helps us to better understand the initial geometry 
and its event-by-event fluctuations. 

In addition, it is important to understand whether 
there is a correlation between different order flow vectors, 
including both the flow planes (direction of the flow vec¬ 
tor) and the harmonics (magnitude of the flow vector), 
and how this correlation between flow vectors modifies 
the underlying p.d.f. One example of the correlation be¬ 
tween different flow planes is the correlation of 4/2 and 
4/3 measured using a 5-particle cumulant: 

2 r,T, 3 _ (cOs(2</3i + 2p 2 + 2^3 - 3(^4 ~ 3yJ 5 )) , , 

u 3l^2} — 3 j (23) 

V 1 

[32] (also denoted by V 23 = v{2, 2,2, — 3, — 3} in [44]). 
Figure 4a shows that uf{ 4 / 2 } from the AMPT model is 
consistent with zero, and since both V 2 and 1)3 are non¬ 
zero, this result proves that there is no correlation be¬ 
tween 4/2 and 4/3 in the AMPT calculations. 

In addition to the correlations between the 2 nd - and 
3 r ’ d -order flow planes we can also study the correlations 
between the other order flow planes. These studies have 
recently gained a lot of attention in the field [20, 42, 44- 
46] . It was proposed to measured these correlations using 
the multi-particle mixed harmonic correlations [ ]: 

(cos(niy>i + ... + n k ipk)) 

= {v ni —v nk cos(ni'I' rai + ... +n k ^ nk )). 

(24) 

These new observables have been measured with V 3 {4 /2 } 
in ALICE [31] and more detailed studies have been pre¬ 
sented in a recent work [42]. Hydrodynamic calcula¬ 
tions [ 20 ] predict that the correlation strength is sensitive 
to both the initial conditions and the details of the ex¬ 
pansion of the system. AMPT simulations, which also 
provide information of initial and final state, will help us 
to understand the role of these correlations in a cascade 
model. 
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FIG. 7: (Color online) Centrality dependence of the final state 
QC{3}(top) and the 2-plane correlations in the initial state 
(bottom) in Au+Au collisions at 200 GeV by AMPT String- 
Melting. 


Using Eq. (24), we can study the n th and m th order 
flow plane correlations in the final state. For example, 
the correlations between ( 4 / 4 ,4/ 2 ) and (4/ 6 ,4/ 3 ) can be 
evaluated via: 

(cos( 4 y>i - 2p 2 - 2^3)) = {vavI cos(44? 4 - 4^2)), 

o (25) 

(cos(6(/?i - 3y?2 - 3 ^)) = (v 6 vl (cos( 64/ 6 - 6T3)). 

As discussed in [ 32 ], these observables can be directly 
calculated in terms of a three-particle cumulant: 

QC{3} 2 n,-n,-n = (cos[n(2yji -ip 2 - ^ 3 )]) 

= [Q2nQ* n Q*n - 2 ■ |Q„| 2 - \Q 2n \ 2 + 2 M\ 
/[M{M — 1)(M — 2)]. 

(26) 


For n =2, we get QC{ 3 } 4 ,_ 2,-2 which is sensitive to the 
correlations between ( 4 / 4 , 4/ 2 ), and analogously for n = 
3 which is sensitive to the correlations between (4/ 6 ,4/ 3 ). 

In AMPT, we can also study correlations between the 
n th - and m th -order symmetry planes, in the initial state 
(4>„, 4> m ). In this model, we can calculate 4>„ assuming 
that the initial spatial energy distribution is proportional 
to the initial spatial parton distribution. Thus, we have: 


$ 


n 


— ATan2 
n 


f (sin {ncj>)) \ 

V(cos (™/o) /’ 


( 27 ) 


where <j> is azimuthal angle of the initial partons. 













For the (<F n , 4> m ) correlation, we can use, similar to 
Eq. (25): 

(cos(4</>! - 2</) 2 - 20 3 )) = (£ A £ 2 2 cos ( 44> 4 - 4 $ 2 )), 

o ( *o) 

(cos(60i - 302 - 30 3 )) = (e 6 eI cos(6$ 6 - 64> 3 )). 

If the initial state symmetry plane <F ra coincides with the 
flow plane 4^, the initial (4> n , <& m ) correlation and the fi¬ 
nal (4^, v F m ) correlation should show a similar centrality 
dependence, and at least should have the same sign. 

We see in Fig. 7 (top) that the ( 4 / 4 , 4/ 2 ) correlation 
has a positive sign, which increases as the centrality in¬ 
creases. In contrast, the ($ 4 , <E> 2 ) correlation is negative 
and decreases with increasing centrality, plotted in Fig. 
4 (bottom). The negative initial ($4, $ 2 ) correlation and 
positive final ( 4 / 4 , 4/ 2 ) correlation observed in the AMPT 
model are in qualitative agreement with viscous hydro- 
dynamic calculations [20] . There is a clear sign change of 
the 4 th -order and 2 nd -order plane correlation during the 
collision system evolution, both in the transport model 
and in the hydrodynamic calculations [20]. On the other 
hand, the 6 th -order and 3 rd -order plane correlation, has a 
negative value in the initial (4> 6 , $3) correlations while it 
is consistent with 0 within uncertainty in the final state 
(4 /q, 4 / 3 ) correlations. In addition, we can see in Fig. 
4 (top) that there is no clear difference of the QC{3} 
from all charged particles (solid symbols) and the like- 
sign (LS) particles (open symbols), which indicates that 
non-flow contributions to the above observables are small 
and the sign change of the plane correlations can not be 
explained by the possible non-flow contributions. 

Considering that the initial anisotropy E n and the final 
anisotropic flow v n are both positive, the above results 
can only be explained if the sign of the genuine correla¬ 
tion, the cosine component, changes during the evolution 
of the produced system. In fact, hydrodynamic calcula¬ 
tions suggest that the final n th -order flow plane 4/„ might 
be not only driven by 4> ra , but might also have contribu¬ 
tions from other symmetry plane(s). The 4^-order har¬ 
monic V 4 and its associated flow plane 4/4 are determined 
by a linear and a quadratic response [47], 

v A e~ li ^ i = w A e~ l4 + u> 4 ( 22 ) e - * 44,2 (29) 

where w A describes the linear response , and W 4 ( 22 ) quan¬ 
tifies the non-linear response. This nonlinear response 
in hydrodynamic calculations couples v A to (V 2) 2 and 
also couples vq to v 2 [ 47 ], and clearly both <F 2 and $4 
contribute to T 4 . Therefore, these results show that 
the initial symmetry plane 4>„ and the flow plane 4f„ 
do not coincident in the AMPT model. Similar re¬ 
sults were also observed in the previous transport model 
calculations using UrQMD model [ 1], and confirmed 
by hydrodynamic calculations [ 49 ]. Thus, the assump¬ 
tion <F n = tFjj used in the so-called true flow calcula¬ 
tions (or symmetry plane flow), ty l {4>„} [16, 50, 51] is 
not valid. In addition, we notice that a new method, 
named scalar product method [ 45 ], is proposed to mea¬ 
sure the flow plane correlations without contributions of 


anisotropic flow. However, it is important to empha¬ 
size that this m ethod is based on the assumption that 
{v 2 V 2 n ) = { v n)V( v 2 n )• f ac t, th e recent ‘standard can¬ 
dle’ calculations SC(rn,n, —m, —n) v [52, 53] have shown 
that there are strong (anti-)correlations between event- 
by-event fluctuation s of v n and v m , therefore the equal¬ 
ity (v 2 V 2 n) — ( v n)V( v 2n) assumed in the scalar product 
method does not hold. Thus, in this paper, we discuss 
the symmetry plane correlations using only the mixed 
harmonic correlations. 
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FIG. 8: (Color online) Centrality dependence of correlations 
between the n th - and m th -order harmonics in the final state 
(top) and initial state (bottom) in Au+Au collisions at 200 
GeV by AMPT StringMelting. 


Not only the correlations between symmetry planes 
but also the correlations between different flow harmon¬ 
ics can be investigated via the multi-particle correlation 
technique. The observable SC(m,n, —m, —n) v was pro¬ 
posed as an unique tool to probe the correlations between 
different orders of flow harmonics and by design indepen¬ 
dent of the symmetry planes [52] . 

In Fig. 8 (top) we see a clear non-zero value for both 
SC(A, 2 , — 4, — 2 ) v (red markers) and SC( 3, 2 , —3, — 2 )„ 
(blue markers) in the final state. The positive 
results of SC( 4, 2 , —4, — 2 )„ and negative results of 
SC( 3, 2 , —3, —2)„ are observed for the presented cen¬ 
trality classes. These (anti-)correlations are more pro¬ 
nounced in peripheral collisions. It indicates that finding 
V 2 larger than V 2 in an event enhances the probability of 
finding V 4 larger than v A and, in addition, the probability 
of finding V 3 smaller than u 3 in that event, as was shown 
in a previous study [ 52] . Also we investigate the correla¬ 
tions of the n th - and m th -order harmonics in the initial 
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state, calculating: 

SC(m, n, -m, -n) e = {e 2 m e 2 n ) - (4J(4> ( 30 ) 

The result is presented in Fig. 8 (bottom). We see a 
positive and increasing trend for SC( 4,2, —4, — 2) e while 
a negative and decreasing trend of 5(7(3,2,—3,—2) E 
is observed. They capture the rough trend of 
5(7(4, 2, —4, —2)„ and 5(7(3, 2,—3,—2) t , in the final 
state but cannot quantitatively describe the central¬ 
ity dependence. This suggests indeed a correlation 
between the initial state SC{rn,n,—m,—n) e and final 
state SC(m,n, —m, —n) v , but this might not be the 
only contribution to the final state. In addition the 
system evolution, might also modify the strength of 
SC(m, n, —m, —n) v . In a previous study [52], three con¬ 
figurations of the AMPT model, named: (a) 3mb; (b) 
lOmb; and (c) lOrnb no rescattering, have been investi¬ 
gated to better understand the SC(m,n,—m,—n) v cal¬ 
culations. In general, the configuration of (a) 3mb gener¬ 
ates weaker partonic interactions during the system evo¬ 
lution compared to (b) lOmb. And the hadronic interac¬ 
tions are turned off for (c) lOmb no rescattering. The 
comparisons between these three configurations shows 
that the value of SC(m,n, —m, —n) v depends on both 
the partonic and the hadronic interactions. These re¬ 
sults suggest that the sign of SC(m, n, —m, —n) v is de¬ 
termined by the initial state SC(m, n, —to, —n) e and the 
magnitude is modified by the multiple interactions dur¬ 
ing the system evolution. 

As we discussed above, the Elliptic-Power function 
gives a better description of the p.d.f. of each single har¬ 
monic. However, it is an open question at the moment 
how the joint underlying p.d.f. including different order 
symmetry planes and harmonics is described and, ad¬ 
ditionally, if these correlations between different order 
symmetry planes and harmonics modify the single har¬ 
monics p.d.f. The investigations presented here begins to 
answer these open questions. Nevertheless, many more 
measurements between different order symmetry planes 
and harmonics are necessary to reasonably constrain a 
joint p.d.f., and ultimately lead to new insights into the 
nature of the fluctuation of the created matter in heavy 


ion collisions. 


V. SUMMARY 

In this paper we presented the calculations of the ini¬ 
tial and final state anisotropies in Au+Au collisions at 
\J snn = 200 GeV using AMPT simulations. It is found 
that the Elliptic-Power function is the only p.d.f. so far 
which describes the event-by-event distributions of the 
eccentricity as well as the triangularity, the anisotropic 
flow from multi-particle cumulants and the relative flow 
fluctuations very well. In addition, the correlations be¬ 
tween different order symmetry planes and harmonics 
have been investigated. A different centrality dependence 
of these correlations in the initial and final state has 
been observed within the same framework of a transport 
model. This result indicates that both the fluctuations in 
the initial geometry and the dynamical evolution of the 
medium in the final state plays an important roles. It is 
currently still unclear how well the underlying joint p.d.f. 
matches the experimental data as these predictions still 
have to be tested at RHIC and at the LHC. The study 
presented in this paper should help us better understand 
the fluctuations of created matter in heavy ion experi¬ 
ments. 
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